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ABSTRACT 


In practice, the measurement of the covariance function (which 
characterizes a time and frequency selective fading communication 
channel) involves envelope detecting the received signals before corre- 
lating. Envelope covariance is insufficient for an accurate evaluation 
of system performance and rather a complex covariance is needed. 
This paper presents a study of these two covariances for Rician sta- 
tistics. The general properties of the Rician distribution are discussed. 
The relationship between complex covariance and covariance of the en- 
velopes of a complex Rician fading communication channel is derived. 
It includes the well known case of Rayleigh fading. For the assumption 
of symmetric spectrum of the narrow-band random process there exists 
a monotonic relationship which allows the evaluation of the complex 
covariance once the envelope covariance is measured. An extension to 
unsymmetric spectrum is given. The results obtained show that the en- 
velope covariance represents a lower bound of the magnitude of the 
complex covariance. 
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RELATIONSHIP BETWEEN COMPLEX- 
AND ENVELOPE-COVARIANCE FOR 
RICIAN FADING COMMUNICATION CHANNELS* 

by 

Erwin Mondr e 

Goddard Space Flight Center 


INTRODUCTION 

The transmission characteristics of a time- and frequency- spread channel can be specified by 
the time-varying transfer function or its Fourier transform (the channel impulse response). Let 
s(t) denote the received narrow-band waveform resulting from the transmission of a sine wave at 
a frequency f + f 0 (f 0 = center frequency) through a fading channel. Then s(t) may be expressed 
by 

s(t) = Re{H(f, t) expj27r(f +f„)t} ’ (1) 

where H(f, t) is the complex envelope of the received signal. Actually H(f , t) is the time-variant 
transfer function of the medium, measured with respect to the center frequency f 0 . Even if this 
function could be measured, using itto characterize a given channel would be somewhat akin to 
using a sample function of a random process to characterize the process. As Gallager (Reference 1) 
has pointed out, these functions contain too much data to be useful without processing. We are in- 
terested only in finding various averages of these quantities. The averages of most interest are the 
mean value, the variance, and the correlation function (or covariance function). 

Let us assume that H(f , t) is a sample function of an ergodic channel which is also stationary 
in the wide sense— the statistical averages of H(f, t) do not change with changes in f and t (Ref- 
erence 2). Then the complex time-frequency correlation function is independent of the location of 
the time interval and frequency band within which the measurement is made: 

R(t, 0) = E^f - y ’ t) H* (f + ’ t + t)| > (2) 

where E indicates an ensemble average and * denotes the complex conjugate. 

♦The research was accomplished while the author held a National Research Council Postdoctoral Resident Research Assoc iateship sup- 
ported by National Aeronautics and Space Administration. 
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MODIFIED CORRELATION FUNCTION 


Bello (Reference 2) presents an experimental technique for the measurement of the complex 
time-frequency correlation function, using independent frequency standards at transmitter and re- 
ceiver site. For this method it is necessary to have frequency standards with relatively good long- 
term stabilities. Unfortunately, such frequency standards in the millimeter-wave region of the 
electro-magnetic spectrum are difficult to use in earth-to-space experiments, because of the weight 
and prime-power constraints on present-day spacecraft. 

An example of such restriction occurs on the ATS-E Millimeter Wave Propagation Experiment 
planned for launch in mid 1969 (Reference 3). In this case the long-term stability is not expected 
to be better than 1 part in 10 6 per year. Therefore a correlation function will be measured using 
only the envelopes (magnitudes of the complex envelopes) of the received signal. With 

H(f, t) = H r (f, t) + jH. (f, t) = V(f, t) exp t)] (3) 

and Equation 1, we find 

s(t) = H r (f, t)cos [277 (f + f 0 )t] - H ; (f, t) Sin[ 2 w(f + f„)t] (4) 


or 


s(t) = V(f, t)cos [^(f + f„)t + 0 ({ , t)] 


(5) 


The output of an envelope detector is now proportional to 


V( f i t ) = I H( f , t )| = [H r 2 (f, t) + H 2 (f, t)] 


1/2 


( 6 ) 


which is exactly the magnitude of the complex channel transfer function. 
The modified correlation function, denoted by R, is defined as 


R (T, fi) = E< 


| H (f ' § ’ t) «(f + ’ t + r 


(7) 


While the modified correlation function is considerably easier to measure, it does not contain in- 
formation about decorrelation of the phase of the transfer function. 
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The following estimation shows that for an ergodie process the modified correlation function 
is always an upper bound of the complex correlation function 

|r(t, n>| < r(t, n> . (8) 

On the assumption that the randomly varying transfer function H( f , t ) is a sample of an ergodie 
process, the correlation function for a certain frequency difference 0 can be calculated by using 
the time average 


R(r) 


1 ini oqp 

T-a> z 1 



H(t)H* 


( t + r) dt . 


0) 


The validity of Equation 8 can be seen from 


l R ( r )l 


1 im 

T — co 


- f 

2 T J. 


2X H( t ) H ( t + r) dt 


1 f 

— 1 im />np I 
T-co * 1 J 


I H( t ) | |H( t + r) | dt = R(r) 


( 10 ) 


Since the received process may have a nonzero mean (this is always the case if the signal en- 
velope is measured), the function really desired is the covariance rather than the correlation func- 
tion. For a complex random varying process, the two-dimensional covariance is defined (Refer- 
ence 4) by 


C(t, H) 


«(f-T ’ l )- E { H (f- T ’ 0} (f+£ ' t t T)-E{H*(f t | ' t + t)} | 


R(r, O) ~ nij m 2 


where m l and m 2 are the complex mean values 


EiH If - o- ’ t 


4 *( 

s{h( 


m 0 = ElH f + t • t 


)} 

>)}• 


According to Equation 7 the modified covariance becomes 


(ID 


(12) 


C(t, n) = R(r, n> - E^f - ■§ ■ t)|J e||h(£ + -§ • t)|j 


(13) 
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From both covariance functions c and c, the coherence time of the channel, T c , and the coherence 
bandwidth, F c , can be calculated. T c is an interval in r over which |c(r, 0)1 is essentially nonzero; 
F c is an interval in O over which |C(0, fi)| is essentially nonzero (Reference 1). In a parallel man- 
ner T c and F c are defined as associated with C(r, O) . But we cannot say that the modified function 
lead to upper- (Reference 5) or lower-bound estimates on the coherence parameters. This follows 
from the fact that a general expression similar to Equation 8 cannot be found for the covariance 
functions. It is necessary to make certain assumptions concerning the statistic of the channel 
transfer function in order to compare these two covariance functions. For a narrow-band Gaussian 
random process with zero mean (i.e., Rayleigh envelope distribution) the relationship between the 
envelope correlation function and the complex correlation function is well known (References 2 
and 6). 

For a channel with known statistic properties, the different average values can be calculated 
theoretically from the first and second probability distributions. Time averages, on the other hand, 
can only be used for experimental measurement techniques. 


CHANNEL WITH RICIAN STATISTICS 

The communication channel between spacecraft and ground station can be considered as a 
multipath channel model with a major strong stable path and additional weak paths. The stable 
path is a result of line-of-sight propagation. The signal then seems to be the sum of a steady sig- 
nal and a fading component, the latter having the usual Gaussian quadrature components. The sig- 
nal envelope resulting from a transmitted continuous sine wave fluctuates according to statistics 
identical with those for the envelope of a sine wave plus additive bandpass Gaussian noise. This 
envelope statistic is usually called Rician fading (Reference 7). 

If we are looking at the channel transfer function for a certain frequency f 0 , the received 
narrow-band waveform is 


s(t) = Re {h( t ) exp ( j 2^ f „ t )\ • 

1 J (14) 

During an interval where the channel can be assumed stationary, the transfer function H(t) is the 
sum of a steady (time-independent) component a + j/3 and a random fluctuation x(t) + jy(t) . The 
sample functions x(t) and y(t) have frequency components only in a narrow band centered on zero 
frequency (References 6 and 8). The random variables x and y, that refer to the possible values of 
x(t) and y(t), respectively, are Gaussian random variables with zero mean and variance 0 2 . 

The transfer function may now be written as 

H(t) = V(t) exp [j0(t)] , (15) 

where the envelope V(t) and the phase 0(t) are slowly varying functions of time. Formally, such 
a representation need not be restricted to narrow-band signals; the concept of an envelope and 
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phase has significance, however, only when the variations of V(t) and 6>(t) are slow compared with 
those of cos 2n f„ t (f 0 is the carrier frequency). The first-order probability density of the envelope 
(magnitude) and phase of the complex transfer function H(t) can be found (References 8 and 10) as 


and 


Pt (V) 



for V > 0 , 


0 , 


for V < 0 , 


(16) 


P x (2) 



r 

/ 1 

1 


I cos e 0 + jFj 

V 2 ’ 

2 ’ 

- z cos 2 0 Q j 


2 0 a ^ exp (- z sin 2 6 ^ ) , for 0 — 0 — 2^ , 


(17) 


otherwise, 


where 


a 2 + /3 2 
Z ‘ ' 2^ 2 

„ (18) 

0 0 - & - arc tan — * 

l 0 is the modified Bessel function of zero order and 1 F 1 the confluent hyper geometric function. 

Figures 1 and 2 are graphs of Pl (V)— or, more strictly speaking, Pl ( v/fi </>)— and Pl (5) 
for different values of z. As expected, the presence of the steady component shifts the average 



Figure 1 — Probabi I ity density distribution of 
envelope for Rician fading signals. 



Figure 2— Probabil ity density distribution of 
phase for Rician fading signals. 
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and the most probable values of the distribution to larger values of the magnitude. For sufficiently 
strong coherent signals, the most probable value of V coincides with the amplitude + / 3 2 of the 

steady component. No signal, z = 0, gives the 



PROBABILITY THAT THE ORDINATE VALUE WILL BE EXCEEDED 

Figure 3— Distribution of the resultant amplitude of 
a constant vector plus a Rayl eigh -distributed vector. 


Rayleigh distribution for the envelope 


exp 


Pi (V) 


V 2 0 a / 


for V > o 


(19) 


J) , for v < o , 

and the phase distribution becomes uniform 


Pi (°) 


* for 0 < 6 < 2tt 


otherwise. 


( 20 ) 


As the signal strength is increased relative to 
the ms fluctuation, the phases are grouped 
progressively closer about the zero-degree 
phase of the carrier, until for no random sig- 
nal the density becomes a delta function at 
9 = arc tan /3/a (<9 Q = o). 


Figure 3, taken from Reference 9, gives 
the distribution of w/j/a 2 + /3 2 for several values 
of z on graph paper so designed that the distribution of the amplitude of a Rayleigh-distributed sig- 
nal (z = 0) expressed in decibels is represented by a straight line with a -45-degree slope. When 
the steady component becomes large, the distribution of the signal level becomes concentrated 
about the value of the steady component (the 0-db line). 


SECOND-ORDER PROBABILITY DENSITIES 

The following sections deal with the determination of auto-correlation function and auto- 
covariance for a Rician fading channel. An extension to the cross-correlation is mentioned under 
’’Covariance for Envelope Detection— Cross-covariance.” 

To derive an expression for the correlation function it is convenient to calculate the second- 
order probability densities. The joint second-order probability density function of real and imagin- 
ary parts of H(t) is, as shown by Middleton (Reference 10): 




(2 7T) 2 |A| 


1/2 


exp 


[- \ X' A” 1 x] , 


( 21 ) 
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where X is the column vector 


H 

~ a. 


x. 

r i 

1 


l 

H.* 

1 i 

-0i 


yi 

H 

- a. 


x 0 

r 2 

2 


2 

H ‘ 2 

-0 2 




whose transpose is the row vector x'. The subscripts 1 and 2 refer to two time instants t t and t 2 . 
The time difference t 2 - 1 1 is denoted by r. The covariance matrix A of the random variables x p 
y j , x 2 , and y 2 is defined by 


1 

0 

r 

x| 

0 

1 

-X 

r 

r 

-X 

1 

0 


r 

0 

1 


( 22 ) 


The determinant of the covariance matrix is denoted by | A | . The mean square of the random vari- 
ables is given by 


- e{xj 2 } = E{y 2 } = f 

J — I 


w( f ) df 


1 , 2 


(23) 


w(f) is the spectral power density, and f is always measured from the carrier frequency f 0 . The 
following statistical averages are functions of r = 1 2 - 1 1 and denoted by 


E { X 1 X 2 } = E { y 1 y 2 } ~ r ( T ) ‘ 



w( f ) cos ( 2n fr) df , 




" E { x 2 y 1 } " M T ) ' ^ 



w( f ) sin ( 2jt fr) df , 


p 2 (r) - r 2 (r) + X 2 (r) < 1 . (24) 

The corresponding joint- distribution densities of the envelope and phase may be found from Equa- 
tion 21 by means of the transformations 


= (v, *<*,»; 

1/2 

1 

«i 

arc tan 

1 i 

H n 


\ 1/2 

to 

arc tan 

H. 

1 2 

) 







Specifically, we have 


P 2 ( V 1' V 2' 0 l’ O 2-' T ) 




where the jacobian J of this transformation is 


3(V V 2' 0 1’ e 2) 


V 1 V 


2 ' 


(26) 


The regions . H, . H f2 , H t <°°j for the variables in a rectangular coordinate system map 

into (OIV,, V 2 <°°) and (ol0j, 0 2 < 2 v) for V. and 0.. Integrating over all phases 0, and 0 2 from 0 
to 2^ gives p 2 (Vj, V 2 ; r). These integrations can be performed with the aid of Equation Bll. 

To simplify considerably the calculation we may assume that the power spectrum w( f ) is sym- 
metrical about the carrier frequency f 0 . In the Millimeter Wave Experiment there is no reason 
that w(f ) should be an odd function of f. Equation 24 shows that for an even function w(f) the term 
\( T ) becomes zero. For the second-order probability density function of the envelopes, (Refer- 
ence 10) gives: 


(V 


) = 


v, v„ 


'P* (l - p 2 ) 


exp 


[- rh -(v^ + v , 2 )/^ 2 (i-P 2 )] 


00 




v r V 2 > 0 , (27) 


where, from Equation 24, p(r) - r(r) , since \(r) = 0. I m (z) is the modified Bessel function of 
order m, and 


e = 2 for m > 1 . 

m 

The power S/N ratio z is already defined in Equation 18. When there is no steady signal component 
(z = 0), Equation 27 reduces to 


v v 

v i v 2 


P 2 (Vj, V 2 ; t) 04 (j - p 2) 


exp 


[-( V 1 2 + V 2 2 )A 2 1 o[ 


pv, 


'p 2 (l 


-p 2 ). 


for V,. V 2 > 0 • 


(28) 


\ 


\ 
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From these second-order probability density functions, the first-order probability density can be 
found by using the general relationship 


lim p - p 2 

r - co z 1 


(29) 


From p(“) = 0 and Equations 27 and 29, we obtain the result already found in Equation 16. 


COVARIANCE OF THE COMPLEX TRANSFER FUNCTION 

The correlation function of the complex channel transfer function may be determined by solv- 
ing the fourfold integral (Reference 6) 


R(r) 


+ 00 

■ Iffl 


dH dH. dH dH. H. 


TT* 

Ho P: 


(w 


"v T J 


(30) 


where p 2 is given in Equation 21. 

A more convenient way to obtain R(t) is the following: From the definition 

R(t) = e|h (t j ) h* (t 2 )} = e{ Hi H 2 *}' (31) 


with 


Hi = + *i + j(/3i +y t ) . 1 = 1 . 2 , 

®{*i} = E{yJ = 0 , 

the correlation function becomes 

R(r) = E{x x x 2 } + E{ yi y 2 } + j[E{x 2 yi }-E{x iy;! }] + (a t + j/3 t ) (a 2 - j/3 2 ) . (32) 

(Note that Equation 32 is valid for auto-correlation as well as for cross-correlation.) Combining 
Equations 24 and 32 gives 

R(t) = 2^ 2 [r(r)- j\(r)] + (a 1 + j£ x ) (a 2 - j£ 2 ) . (33) 
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With the mean values of the transfer function 


m i " E { H i} ” a i + i^i - 

m 2 = E { H 2 } = a 2 + ^2 ’ (34) 

the covariance, defined as 

C(t) = R(r) - m 1 m* , (35) 

becomes 

C (r) = 2/< 2 [r(r) - j\(r)] = 2 *p 2 p(r) e-’**- 7 "* . (36) 

This result shows that the covariance does not depend on the coherent part of the transfer function. 

The value of C(r) for r = 0 is the variance of the random process. Equation 36 gives, since 
r(0) = l and>t(0) = 0 (compare Equation 24), 


C( 0 ) = 2 >/< 2 . 

The normalized covariance or correlation coefficient becomes 


— - — = t(t) - j\(r) = p(r) 

C( 0 ) 


(37) 


(38) 


where p(t) is the magnitude and 4>(r) the phase: 


<£(t) 


arc tan 



The absolute value of the normalized covariance from which the channel parameters can be de- 
termined is 


p(r) - [r 2 (t) + \ 2 (r)] 1/2 , 0.1 p(r) < 1 . 


(39) 


The complex normalized covariance depends only on the random part of the complex channel func- 
tion and therefore yields equal results for Rayleigh and Rician fading channels. On the other hand 
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(as will appear later) for envelope detection, the correlation function and mean value both depend 
on the ratio of the coherent and random parts of the process. 


COVARIANCE FOR ENVELOPE DETECTION 


Auto-covariance 


The correlation function of the envelope R(t) can be calculated by taking the average value of 
the second-order probability density function, Equation 27. Since this equation is only valid for 
X(r) ■ 0, this restriction is held throughout this discussion of covariance. The integral 


R(t) 


CO 

h 


V 2 P 2 ( V 1- V 2 ; r) dVj dV 2 


(40) 


can be evaluated in different ways, as shown by Middleton (Reference 11). Because his results are 
incorrect, Appendix A gives one method in detail. 

The simplest way follows after expanding i m [pVj v 2 /V (l -p 2 )] according to Appendix B, Equa- 
tion B9, and using term-wise integration with the help of Equation BIO. This yields 


R(-r) 


2P 1 


CO 

(l-p 2 ) 2 e‘ 2 ^ /1+ ^ 2 ^ 


( m !) 2 


/ 1 - p \ m S~' p 2n r 2 (m +n + 1.5) 
' 1 +p ' Z_I n! (m + n) ! 

n = 0 


• ^ (“ n - 0.5; m + 1; -z • (41) 

where 1 F 1 (a; b; z) is the confluent hyper geometric function, defined by 


1 F 1 (a; b; z) - 


az a(a + 1) z 2 

1 + W) + b(bTT) 2i + 


(42) 


The mean value of the received signal envelope may be obtained from the correlation function 
Equation 41, by setting r = oo. With p(°°) = 0, we find 


R(°°) = 2sp 2 r 2 (1.5) X F X 2 (“0.5; l;-z) = m 2 . 


(43) 
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The mean and mean-square values can also be calculated in a more direct way from the first-order 
probability density function, Equation 16. Equation BIO then gives 


E{V> = m 


I V Pl (V)dV = </< |/^ jF, (-0.5; 1; - z) 

Jo 


(44a) 


e{v 2 } = V 2 Pl (V)dV = 2'p 2 ( 1 + z) • (44b) 

Jq 


Equations 43 and 44a yield the same result for m, because f(l.5) - f^/ 2. 

The normalized covariance for the envelopes of the received complex function becomes, using 
Equation 13, 


p( r ) 


C(T) 

C(0) 


R(r) 


(45) 


where 


C(0) = R(0) ~ m 2 = E{v 2 } - E 2 {V} = a 2 . 

The correlation function, the mean value, and variance are defined in Equations 41, 44a, and 44b, 
respectively. Equation 45 shows a relationship between the normalized covariance 3 of the enve- 
lopes of the complex process and the magnitude of the normalized complex covariance p (defined 
in Equation 39). 

Since this relationship is rather complicated, we develop Equation 41 in a power series of p 
up to p 2 to obtain an approximate value. This is also to compare the results with those calculated 
by Lawson and Uhlenbeck (Reference 6), where the correlation function is given also as a power 
series up to p 2 . Combining 


R - A Q + Aj p + A 2 p 2 + 


with 


A; = IT 


1 d 1 R 


dp 1 


'p~0 


i r 0, 1, 2, ’ * * 
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and Equation B3, gives 

A 0 = m 2 (compare Equation 43) 

A, = z[| jFj (- 0.5; 2; -z) - jFj (-0.5; l; - z)] 

A 2 ,= {[| jFj (- 0.5; 2; -z)] - § ( 1 + z) 1 F 1 (- 0. 5; 2; - z) jFj (0. 5; 2; - z) 

+ \ (z 2 +2z + (0.5; 2; - z)} • 


The result is still too complicated; it is convenient to consider separately the cases where z < 1 
and z > 1. With the developments for the confluent hypergeometric functions, Equations B1 and B2, 
for small coherent parts of the transfer function (z < 1) Equation 45 becomes: 


p 



< 


l, 


P < l • 


For large coherent parts (z > 1) Equation 45 becomes 



1 4z ~ 8z 2 


z > 1, p < 1 • 


(46) 


(47) 


In Reference 6, the approximations for the correlation function show only the first term in each ex- 
pression proportional to p and p 2 , respectively. Figure 4 plots p versus p, as in Equations 46 and 
47, for different values of z. The curves are valid only for p < 1, because only terms up to p 2 are 
taken into account. Due to the normalization, p must become 1 for p = 1 (which corresponds to 
r = 0) and for all values of z. For the covariance, it is always true that 


P < P , (48) 

whereas, for the correlation functions, R — I r| (see Equation 8). 

For an exact calculation of r (t) Equation 41 is not very useful because it converges less as p 
increases. Since the term (l - p 2 ) becomes zero for p = 1 (r =0), the sum must become infinite 
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Figure 4— Normalized covariance of envelope versus 
normalized complex covariance for Rician fading. Ap- 
proximation valid for p < 1. 


to provide a finite value for R(0). For this 
reason, another way of evaluating the inte- 
gral in Equation 40 will be used. This is 
achieved with the aid of the transformation 
(Reference 11) 

Vj = A yi ~ P 2 /w e ^ /2 

V 2 - ^ )/l - p 2 /w e~ M/2 , (49) 

followed by the successive application of the ex- 
pression for the product of two Bessel functions 
and using the integral form of the modified 
Bessel function of the second kind. Appendix A 
evaluates this integral. The correlation func- 
tion in this case becomes 


OD 00 


R(r) = 2P 2 e ~ 2zn 'P L L 


p m z m + n 


m - 0 n= 0 


l! (1 + p) 2m + 2n 


n 

L 

k=0 


T(m +n -k + 1.5) T(m+ k + 1.5) / . . 

" ..F.l-k-O.S, - n + k-0.5;m + l;p 2 ), f50'l 

k!(n-k)!(m+k)!(m+n~k)! 21 ' ' \ v) 


where 2 f % (a, b; c; z) is the hypergeometric function, defined by 


V"> ( a )k ( b )k b 

2 F l ( a ’ b; c; z > = /_> _ (c) k k! z • 


with 


(a) k = a(a + l) (a+2)-“ (a+k- 1) , and (a) Q = 

For the limiting case of zero steady component (z = 0), Equation 50 reduces to 
R 0 (r) a 2^T 2 (1.5) 2 Fj (-0.5, -0.5; 1 ; p 2 ) . 


(51) 


(52) 
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The hypergeometric function 2 F t may now be expressed in terms of the complete elliptic integrals 
of the first (k) and second (e) kind: 


2 F, (-0.5, -0.5; 1; p 2 ) 


77“ EC/ 3 ) 77 (l P 2 ) K(P) 


(53) 


We find 


R 0 (t) = 02 \ jz {p) _ (1 .^2) K(p) ] f (54) 

a well known relationship for Rayleigh- envelope fading characteristics (References 2, 6, and 12). 

In this case the mean value, by Equation 44a, becomes 

m 0 = 0 /V2 , (55) 


and the variance, by Equation 44b, becomes 

a 0 2 = e{v 2 } - in 0 J = 0 2 ( 2 “ ^ ) (56) 

These values can, of course, also be achieved from the general formulas: 

R(c°) - m 2 

R(0) = E{V 2 } • (57) 

From p(0) = 1, pC 03 ) = 0, and Equations B7 and B8, Equation 52 yields 

r o (°°) = 0 2 n/2 

R 0 (0) = 20 2 . (58) 

Figure 5 plots R (that is, r/ 20 2 ) versus p for different values of z. The normalized covariance 
p for the Rician envelope statistic as a function of the magnitude of the complex normalized co- 
variance p y as defined in Equation 45, is calculated from Equations 50 and 44a and b for different 
values of z and shown in Figure 6. For Rayleigh distribution (z = 0), compare Schwartz (Refer- 
ence 7, p. 479). 

As expected, p and p differ most for the Rayleigh statistic (z = 0). This difference decreases 
with increasing coherent part of the received wave form (increasing z). For z -oo ? we find p = p. 
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Figure 5— Correlation function of 
envelope versus p(r). 



Figure 6 — Normalized covariance of envelope versus nor- 
malized complex covariance for Rician envelope statistics. 

According to the definition of p in Equation 45, 
this value must become 1 for p = 1, which ob- 
tains when delay r = 0, since a- 2 = R(0) -m 2 . Be- 
cause of the monotonic relationship between p 
and p it is possible to determine p, once p and 
the parameter z are measured. 


The parameter z can be calculated by measuring the first moment (i.e., the mean value or dc 
component) and the second central moment (i.e., the mean square of the varying component or ac 
component) of the recieved signal envelope. Equations 44a and b show that both are functions of \p 
and z: 


m = E{V} = 4>-fp2 jFj (- 0.5; 1; - z) 

& 2 = E { v 2 } - m 2 = </- 2 { 2(1 + z)- \ jFj 2 (- 0 . 5 ; 1 ; - z)} • ( 59 ) 


Once the ratio m 2 / 2 " 2 (which is only a function of z) is calculated, the parameter z can be determined 
with the aid of the following equation: 

m 2 _ f 4(1 + z ) 

2& 2 " 2 L TT jF 2 (-0.5; 1; - z) 

Figure 7 is a graph of Equation 60. 
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In case of Rayleigh fading (z = 0), 


m 2 0.5 



and, for z -a>, Equation 60 becomes 

2 

m 

lim - z + 0. 75 • 

' z -CO 2cr^ 

Since a- 2 is the mean noise power at the 
output of the detector andm 2 /2 the mean signal 
power (1-ohm resistor), m 2 /2o- 2 is simply the 
power S/N ratio at the output of the envelope 
detector. On the other hand, the parameter z 
is a characteristic of the channel transfer 
function. It is defined (compare Equation 18) 
as the ratio of the mean-square values of the 
steady and fading components. 



Figure 7— Output-power S/N ratio versus input-power 
S/N ratio for Rician fading and envelope detection. 


Cross-covariance 

Here the results of the previous sections are extended to the correlation between two signals 
at different frequencies. Now we must calculate the cross-correlation function of two signals that 
do not necessarily have equal coherent signal components. Since the normalized complex covari- 
ance does not depend on the coherent components, Equation 38 is still valid for the correlation 
between two different signals. The modified correlation function (correlation function of the enve- 
lopes), on the other hand, depends on the coherent parts. This is because the second-order prob- 
ability density for the envelopes is a function of the two parameters and z 2 — the s/N ratios of 
the two different received signals (Reference 10), and is given by 


P 2 (V,,V 2 ; t) 


V ! V 2 

(l~P 2 ) 


exp 


V 2 + V 2 
v 1 v 2 


20 s (l-p 2 ) 1 - P 


Z l + Z 2 r— 

+ -7T7 


ir T ( pv * v 

Z_, (i -/ 


Y2 

I "> \W 1 


MMt 


2 1 -p 2 j 

i2 fc-pfc 


i- p 2 


(61) 


To obtain the modified cross-correlation function, we can use the method outlined above for 
the auto- correlation function. The integral in Equation 40 in conjunction with Equation 61 can be 
solved in the way that is shown in Appendix A. It is convenient to use the two parameters 


Z 0 = V z i Z 2 
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y 


(62) 


/V* 


instead of z x and z 2 . The cross-correlation for the envelopes becomes 


R c ( T ) - 2</> 2 exp 


n 00 oo 


y + l/y - 2py r 

1-P 2 J Z__. Z-. £m ra! z ° 


m=0 n= 0 


[( 1 “ p/y) ( 1 - py)] m [(7 “ p) 2 /y] / I - py \ 2k T(m +n - k + l.5)r(m+k+i.5) 

( 1 - ^ 2 ) 2m+2n L_ i ' ^ ~ P / k ! (n-k)! (m + k) ! (m + n - k) ! 

v ' k = o 


* 2 F 1 (-k - 0.5, ~n +k -0.5; m+ 1; p 2 ) . (63) 

We obtain the relationship between the normalized complex cross-covariance p and the normalized 
cross-covariance of the envelopes p c (r) by 


R c (r)-m 1 rh 2 

M T > = ' ara 2 ’ (64) 

where R c (r) is given in Equation 63, and the mean values m. and standard deviations 6-. are speci- 
fied in Equation 44. The subscript values i = 1, 2 refer to z x and z 2 , respectively. 



Figure 8 — p(r) versus p(r) for z Q = 0.5 
and different values of y. 


It should be mentioned that y = 1 repre- 
sents the case z 1 = z 2 = z Q (e.g. auto-correlation 
The parameters y and z Q may be calculated 
from Equation 62 by measuring the dc and ac 
components of the two received signals and us- 
ing Figure 7. Figure 8 is a graph of Equa- 
tion 64 for z Q = 0.5 and different values of y. 

As expected, the curves for given values of y 
and l/y are exactly the same. This can also 
be shown by replacing y by l/y in Equation 63. 
This leads to 

R c (t, y) s R c (t, l/y) . 

The modified covariance is highest for y = 1 
(auto-covariance) and becomes smaller if y 
decreases or increases (different steady com- 
ponents of the two received signals). Note that 
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cross-covariance is normalized with respect to the variance of the two received signals and not to 
some carrier variance. 

Channel Parameters 

The preceding consideration shows that for Rician fading channels the following expression is 
always valid (compare Figures 6 and 8 and see Equation 48): 


P ( T ) - P.( T ) 


(65) 


The difference p-p becomes smaller with increasing coherent signal component (increasing z). 


If the modified covariance functions are used to calculate the coherence parameters of the fad 


ing channel (i.e., no correction according to the 
tions mentioned above is used)j the coherence 
time t and the coherence bandwidth F n are al- 
ways less than the values calculated from the 
complex covariance: 

T c < T c 

F=iF t ' (66) 

Figure 9 is a 2-dimensional section 0°. T ) 
of a 3-dimensional figure with axes: p, t, and 
a (it shows the case for a = 0). T c and T c are 
shown, but not F c and F c (which would require 
a p, n graph). 


correspondence between the two covariance func- 



Figure 9 — Definition of coherence time and relationship 
between modified and complex covariance. 


On the other hand, the frequency spread B and the time spread L of the channel are wider than 
their actual values 


B > b , 


L > L . (67) 

since they are derived from the scattering function, which is a Fourier transform of the two- 
dimensional normalized covariance. 

Few experiments seem to have been made to verify the usefulness of the theoretically derived 
channel parameters, particularly in the millimeter region of the electromagnetic spectrum. If the 
results achieved from the MillimeterWave Propagation Experiment on board the spacecraft ATS-E are 
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corrected as outlined in the previous chapters, more useful data about the time- and frequency- 
selective fading properties of the channel between satellite and ground terminal are obtained. Al- 
though the stability of the instrumented oscillator allows only the measurement of the modified 
covariance, the magnitude of the complex covariance for the channel-transfer function can be cal- 
culated if the received random signal has the expected Rician probability distribution. From the 
magnitude of the two-dimensional complex covariance function, the coherence parameters of the 
channel are characterized and the communication- signal propagation effects of the channel can be 
determined. The corrected data represent a better description of the fading channel than the data 
achieved from the measured modified covariance, and allow better and more accurate planning and 
design of future communication links through the atmosphere for millimeter waves. 


RANDOM SIGNAL WITH UNSYMMETRIC SPECTRUM 

The modified correlation function in the preceding sections was calculated for a symmetric 
power spectrum. In the more general situation where A ,(r) ^ 0 (spectrum not symmetrical in the 
narrow band about center frequency) the second-order probability density function is much more 
complicated. By introducing polar coordinates in Equation 21 and integrating over all phases (with 
the aid of Appendix B, Equation Bll), we finally obtain 


P 2 {V 1 , V 2 ; t) 


^1 ^2 2z( 1 - p cos <£) 

r (1 -r>j exp L ■ 


V l 2 + V 2 2 

2 0 2 (l -P 2 ). 


^ ( PV 1 V 2 \ 

• *■ ( ^— ft - P 2 ) ) I. ( 1 W, ) I, (bV 2 ) cos (#-»), for v l ,v 2 >0, 


( 68 ) 


where 


and, from Equation 24, 


z - (a 2 +/3 2 )/20 2 , 


+ p 2 ~ 2p cos 0 
0 1- p 2 


€ 


arc tan 


p sirup 
p cos <£ - 1 


P( T ) ‘ [r 2 (r) + X 2 (r)] 1/2 

X(r) 

<£(t) - arc tan • 


(69) 
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Equation 68 can be found in Middleton (Reference 10) but there is a sign error that leads to a cos0 
term instead of cos (2£-<£) . For \ = 0, Equation 68 reduces, of course, to Equation 27. 

To obtain the auto-correlation function of the envelopes, we use the same procedure as shown 
in Appendix A for the case A. = 0. This leads to 


K ( r ) 2 ^ 2 


exp 




2 z 


1-/0 cos0 ~l Y"* y- p m z m + n (l + p 2 - 2p COS <fi) n 


1 -P z 


m ~ 0 n=0 


m! (l “ p 2 ) 


2m + 2n 


* COS 


n 

[m(2 €~<t>)] ^ 


r(m + ri“k + 1 . 5 )r(m+k +1.5) 
k! ( n - k ) ! ( m + n - k ) ! ( m + k ) ! 


2 F 1 (~k -0.5, ~n + k -0.5; m+ 1; p 2 ) 


(70) 


The correlation function for a symmetric spectrum is a special case of Equation 70, and we ob- 
tain Equation 50 by setting 0 = 0. Since the first and second moments do not depend on 0, they are 
the same as shown in Equations 44a and b. 

The graphs showing the normalized modified covariance versus complex covariance are similar 
to those in Figure 8. The modified correlation function and therefore the covariance of the enve- 
lopes, becomes smaller with increasing <t> (or A.); this is due to the decorrelation effect of the asym- 
metric spectrum. Since is an even function of 0 , the same results are obtained for ±0 or dA. 

A numerical calculation shows that the difference between the curves for 0 = 0 and |0| - tt/S is 
less than 7.5 percent for all values of p. 

If the spectrum is asymmetric about the center frequency, the complex covariance cannot be 
calculated by measuring of p(r) alone, since the angle as well as the absolute value of the complex 
covariance are unknown. The envelope detection yields insufficient information about the complex 
correlation function. Only for the Rayleigh fading channel (z = 0), the above equations become in- 
dependent of 0 . Therefore, in this case there is only one relationship between p and p valid for all 
values of 0. In Rician fading channels (z ^ 0) we must assume that \ = 0 in order to have a definite 
relationship between p and p. As already mentioned, this assumption is quite reasonable for the 
received narrow-band signal. 


CONCLUSIONS 

This paper describes the properties of the modified correlation functions and presents a com- 
parison of the modified and complex correlation function for the Rician fading channel. The modi- 
fied correlation function, i.e., the correlation function of the envelope of the received signal, is 
considerably easier to measure than the complex function, and, in cases where the long-term sta- 
bility of the measurement system oscillators is not good enough to determine the complex cor- 
relation function, is the only function that can be measured. 
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It has been shown that for Rician fading statistics the modified covariance is a lower bound of 
the magnitude of the complex covariance. This leads to a useful estimate of the coherence param- 
eters for the fading channel, as determined by the modified functions (under "Channel Parameter"). 

Because the relationship between p{r) and p(r) is monotonic and single-valued if the power 
spectrum of the received narrow-band process is symmetric about the center frequency, we can 
determine the magnitude of the complex covariance p(r) by measuring the modified covariance p(t). 
To test the equivalence of the probability density function for the sampled data to the theoretical 
Rician probability density function the so-called chi-square goodness-of-fit test described by Bendat 
and Piersol (Reference 13) can be used. As pointed out under "Auto- covariance," the s/N ratio z 
can be calculated, once the mean value and the mean-square value of the received signal are meas- 
ured. Since, in the relationship between p(r) and p(r) (Equations 45 and 50), p(r) is the independent 
variable, a numerical approximation method or the graphs in Figure 6 can be used to calculate p(r) 
from the measured value p(r). As expected, the difference between p and p becomes smaller with 
increasing coherent part of the signal. The greatest difference can be found for zero coherent part, 

i.e., Rayleigh fading. The coherence parameters can then be found from the magnitude of the com- 
plex covariance p(t), and provide a better knowledge about the properties of the fading channel. 

For design of future communication links, these corrected data are more useful than those calcu- 
lated directly from the modified covariance. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, May 15, 1968 
160-21-10-01-51 
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Appendix A 


Evaluation of Integral of Equation 40 


To obtain the auto-correlation function for the envelopes of the received random signal, we 
must solve the integral 


R(r) 


00 

JJ 


V 1 V 2 1 


(V 


V 2 ; r) dVj dV 2 


where p 2 is given in Equation 27. To evaluate this integral we use the transformation 


(Al) 


v, = 


2 w l /2 - m /2 


'pyi - p 


2 w 1/2 e -M/2 


(A2) 


where 0 < w <co and - a> < M < oo. The jacobian of this transformation from V. and to w and p is 


*(v y 2 ) 

5(w, M) 


= \ '/> 2 (l ~ P 2 ) 


(A3) 


In conjunction with Equation A2, the arguments of the Bessel functions in Equation 27 become 


^ V 1 V 2 
0 2 (l -P 2 ) 

Y57 v x 

</- (i + p) 

/2I v 2 

4 K 1 + p) 


pw , 


sw 1//2 e M/2 


sw 1/2 e ~M/2 
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Zl-P 2 
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(A4) 
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In addition, we need 


v x V 2 = (I -p 2 ) W , 

v 2 2 + v 2 2 

2^ 2 (i- p 2 ) = WCOsh ^ • ( A5 ) 

The integral in Equation A1 therefore becomes 


R = -j \p 2 (l - p 2 ) 2 e" 2z/1+p ^ el dw I d^w 2 e -wcosh ^ * I m (pw) I m ( sw 172 e^ /2 ) I m ( sw 1/2 e _ ^ /2 ) *(A6) 

J 0 J-CD 


The product of two modified Bessel functions can be written with the aid of Watson (Reference 14, 
p. 148, Equation 2) as 


0 


i; 1/2 e M/2 


KO 


h 2 V co /i 

. W w / V i l 9 sv 


1/2 e ^/2 


(m + n)! 


2 Fj (-n, -m -n; m+ 1; e" 2 ^) . ( A7 ) 


For the hyp ergeo metric function we use the power series representation, Equation 51: 


2 F 1 (""■ 


00 



k= 0 


(~ n ) k (~ m ~ n ) k 
k! (m + l) k 


e' 2k m . 


(A8) 


Since n is an integer, (-n) k becomes zero for k >n + 1 . Therefore the series in Equation A8 con- 
sists of a finite number of terms. The upper limit in the sum can be replaced by n. 

The integration over the variable fi, which may be written in the following form, if we use 
Watson (Reference 14, p. 182, Equation 7) yields 


n^-wcosh^ = (w) , (Ag) 

where K 2k _ n is the modified Bessel function of the second kind and order ( 2 k -n). Now we must 
evaluate the integral over w. Using Watson (Reference 14, p. 410, Equation 1) we find 



1 “ 

•'0 


dww 2+m+n I m (pw)K 2k _ n (w) = p' 


r(m + n ~ k + 1 . 5) F(m + k + 1 . 5) / , , o\ 

2 -i-.-n r(nl+1) • 2 F 1 (m + n-k+1.5,m+k + 1.5;m+l;p 2 ). (A10) 
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With the relation 


2 Fj (a, /?; y; x) = ( 1 - x)^ a P 2 F t (y - 0, y ~ a; y; x) , 

and, by suitable rearrangement of the above formulas, the auto-correlation function becomes 
finally 


00 CD 

» r-i p m z m+n 


m! (1 +/5) 2m+2n 


y r ( m + n - k + 1 . 5 ) r ( m + k + 1 . 5) { _ k _ 05> - n + U - 0 . 5; m + 1; ) 

/ . k!(n-k)!(m+k)!(m+n-k)! 21 


Note that in Reference 11 the factor m! as denominator is missing. 


(All) 


(A12) 
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Appendix B 


Special Functions and Integrals 


This appendix is an assembly of some mathematical results used throughout this report. Other 
functions and integrals, already mentioned, will be omitted. 

Hyper geometric Functions: 

The confluent hypergeometric function is defined by 


(a) n * n 

iF, b; z) “ (b)n n! ’ (Bl) 


where (a) n = a(a + l) (a+ 2) • " (a + n - i), and (a) 0 = 1. The asymptotic series for l F 1 (a; b; z) for 
large values of z is 


l F i (a; b; - z ) 


y-n ( a ) m ( a b + 1 )n 
/ I ml z m 


F (b) z ~ a 

r(b-a) ^ m! z 

m= 0 


r(b)z 

r<b 


- a) L 


a( a - b + 1 ) 


(B2) 


There exists a set of recurrence relations, which have been used in the major section on covariance: 


- ZF 11 

= bF oo ■ bF 10 - 


aF n 

= (a-b)F 01 + b F 00 , 


abF ,o 

= b( a +z)F 00 - z(b -a)F 01 , 


( b ~ a) zF 01 

= b(z +b-l)F 00 + b(l-b)F 0 _ 1 , 

(B3) 


where the symbol F k ,^ is an abbreviation defined by 

F u = i F i ( a + k ; b + *'< z ) • 
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The relation 


e z jFj (a; b; z) - ^ (b- a; b; -z) 


is also needed in the discussion of covariance. 

The Gaussian hypergeometric function is represented by 


V"* ( a ) n ( b ) n 

2 F i < a ' b; z) = ^ "(c) n ri! ■ ■ 


1*1 * l 


Useful relations are 


2 Fj (a, b; c; z) - ( 1 - z ) c a b 2 Fj ( c - b, c - a; c; z ) 


2 Fj (a, b; c; 1) 


T(c) T(c - a^Jb) 
f(c - a) P(c - b) 


Re (c) / 0, -1, “2, 


Re{c - a - b} / 0, “1,-2, 


2 Fj (a, b; c; 0) = 


Modified Bessel Function: 


The modified Bessel function (Bessel function of purely imaginary argument) of the first kind 
is defined by the power series 


* viiiP. 

( z ) / n! (m + n)! 


The following integral, used in the discussion of covariance, can be derived from Hankel’s ex- 
ponential integral; compare Watson (Reference 14, p. 394). 


: iV [ exp (£)] iFi ( 2 
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To solve an integral of the form 



(a cos 0) dd , 


we can use the following series development: 


exp (a cos 6) 


L 

m= n 


£ m I n. ( a ) cos ( m 0 ) ■ 


(Bll) 
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